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Abstract 

We show that there exist maximal globally hyperbolic solutions of the Einstein-dust 
equations which admit a constant mean curvature Cauchy surface, but are not covered by 
a constant mean curvature foliation. 
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1. Introduction 

The most widely advocated candidate to be used as a choice of time in cosmo logical 
spacetimes, i.e., spacetimes which are 

a) globally hyperbolic with a compact Cauchy surface, and 

b) solutions of Einstein's equations coupled to 'reasonable' (see Sect. 2) non-gravitational 
fields 

is that which is specified by a constant mean curvature ('CMC') foliation. Supporting this 
advocacy is the belief that a generic cosmological spacetime admits a global CMC foliation 
which is unique up to leaf-preserving diffeomorphism. While this belief may turn out to 
be true, it was shown some years ago [1] that there are some cosmological spacetimes 
which do not admit any CMC Cauchy surfaces. Here we show further that there are also 
cosmological spacetimes which admit CMC Cauchy surfaces but do not admit a global 
foliation by such hypersurfaces. 

The spacetimes which we show have this property are quite special *: the matter 
fields in them are pressureless fluids, i.e. dust, which are contrived to form 'shell-crossing 
singularities', and the spacetime fields are all invariant under a T 2 -isometry group. While 
the latter feature is merely a convenience, it may well be that the formation of a shell- 
crossing singularity is crucial to the failure of the CMC slicing to cover the spacetime. 
This issue needs further study. In any case, we believe it is useful to know that there are 
cosmological spacetimes in which CMC-based time exists for a while, but does not cover 
the whole globally hyperbolic spacetime. 

Our results here, to an extent, build on previous work by one of us[12]. In that work, 
it was shown that there are cosmological spacetimes (with T 2 isometry and satisfying the 
Einstein-dust equations) which, starting from initial data with finite constant mean cur- 
vature To, develop a (shell-crossing) singularity arbitrarily soon - To + e - in CMC time, 
provided the solution exists that long. What we show here is that, while the singularity 
stops the CMC foliation from continuing past tq + e, for some (generic?) choices of ini- 
tial data, the spacetime continues to develop globally hyperbolically in regions which are 
spatially distant from the singularity. Hence, there are spacetime regions which are left 
uncovered by the CMC slicing. Note that while examples of spacetimes with incomplete 
CMC foliations can be constructed trivially by cutting pieces out of certain spacetimes 
(see, for example, p. 472 of [13]), such examples are not maximal globally hyperbolic 
developments of initial data, as are our examples described here. 

To prove the existence of spacetimes with these properties, we first (Sect. 2) briefly 
discuss cosmological spacetimes satisfying the Einstein-dust equations, noting especially 
results concerned with stability, and we describe the special form of the fields and some of 
the field equations if we assume that the spacetimes are T 2 symmetric. Next, in Sect. 3, 
we describe the special sets of initial data which were used in [12] to produce singularities 
in arbitrarily short CMC time. Then in Sect. 4, we specify the 'quiet zone' restrictions on 
the data which we use to produce our present results, and we state and prove the results 
in Theorem 3 in that section. 

* Note that, thus far, the spacetimes which are known not to admit any CMC hyper- 
surfaces are also special in that they are solutions of the Einstein-dust field equations. No 
vacuum solutions are known to have this property 
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2. Cosmological spacetimes with dust 

As noted earlier, a spacetime (M 4 , g a p-> ^) with metric g a p and non-gravitational fields 
\& is called cosmological if (M 4 , g) is globally hyperbolic with compact Cauchy surface, and 
if g a p and \1/ together satisfy the coupled Einstein-matter system: 

G af3 [g] = 87rT al3 [g,y} (la) 
.F[<7,*] = (16) 

where T a p[g, ty] is the stress-energy tensor of \1/ and g a /3, and where T\g, \P] = represents 
the field equations which the chosen field theory imposes on g a p and The globally 
hyperbolic condition implies that M 4 = E 3 x R for some compact three-dimensional man- 
ifold E 3 , and also implies that any spacelike embedding of E 3 in (M , g) is a Cauchy sur- 
face[2]. Global hyperbolicity does not, however, guarantee that the spacetime (M , g a p, \&) 
is uniquely determined by initial data on a Cauchy surface. For this property to hold, the 
PDE system (1) must be hyperbolic in an appropriate sense [6]. This PDE hyperbolicity 
condition, along with an energy condition (weak, strong, or dominant) [7]) on the energy- 
momentum tensor, is what we mean by 'reasonable' non-gravitational fields in the definition 
of cosmological spacetimes. 

Our focus in this paper is the Einstein-dust field theory, in which the non-gravitational 
fields consist of a non-negative function /i (the proper energy density of the matter) and 
a unit vector field u a (the four- velocity of the matter). The stress-energy tensor for this 
theory takes the form 

T Q/3 = (lU a Up, (2) 

and the non-gravitational field equations are just those derived from the Bianchi identity 

= V a G a p [g] 
= V a T a p\g,fji,u] 

One readily verifies from (2) that the weak, strong and dominant energy conditions are 
all satisfied by the Einstein-dust theory. In contrast, it is not so easy to verify that the 
PDE system for this theory is hyperbolic. The procedure of reduction to a symmetric 
hyperbolic system, which straightforwardly shows hyperbolicity for the Einstein equations 
coupled to many matter fields, runs into difficulties in the case of dust. (In contrast it 
does work for the Einstein-perfect fluid equations with a reasonable equation of state and 
strictly positive pressure.) By rewriting (l)-(3) in the equivalent form 

-u 7 V 7 G Q/ 3 = - / u(V 7 -u 7 )'u a 'u / 3 (4a) 
-u 7 V 7 w/3 = (46) 
V 7 (^ 7 ) = (4c) 

Choquet-Bruhat is able to show that the Einstein-dust field equations are hyperbolic in 
the Leray sense [3], [9]. 

It is important for our present purposes to know not only that one can determine the 
spacetime (M 4 , g a p, /i, u a ) from knowledge of initial data on a Cauchy surface E 3 o , but also 
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that at least for a finite (in time) spacetime neighbourhood M 4 tq _ a To+X ) °f £v > the fields 
in this determined spacetime (M^ tq _ x To+x yg a p, ii,u a ) are stable under perturbations of 
the initial data. The theorem which guarantees this stability for PDE systems which are 
hyperbolic in the sense of Leray was proven by Choquet-Bruhat in [4] using harmonic 
coordinates. Applied to the Einstein-dust equations (4) this result says the following. 

Proposition 1. ( Choquet-Bruhat [4]) Let (Mq , g^g, Ho, Uq) be a cosmological spacetime 
satisfying the harmonically reduced Einstein-dust equations and let (Mq, g^pi Ao> u q) be 
restriction of this solution to a region Mq bounded between a pair of Cauchy surfaces in M 4 . 
Let (g^g, fiQ, Uq) be contained in the Sobolev space Hi(M 4 ), and for a fixed Cauchy surface 

E 3 in M 4 , let the corresponding Cauchy data Co be contained in ^(E 3 ). There exists an 
open set W C ^(E 3 ) with Cq G W such that every set of initial data C G W generates a 
cosmological spacetime (M 4 ,g a p, fa,u a ) which satisfies the harmonically reduced Einstein- 
dust equations, and also (g a f3, fa,u a ) Hi(M 4 ) . Moreover the mapping which takes a point 
ofW to the corresponding solution is continuous (in fact differentiate) with respect to the 
Hj topology. 

The spacetimes that we show here have an incomplete CMC slicing are Einstein-dust 
cosmological spacetimes, and in addition have E 3 = T 3 and admit a T 2 isometry group 
which acts spatially. In a neighbourhood of some fixed CMC initial slice Ejj, the metric 
for such a spacetime can be written in the formfll] 

g = - a 2 dt 2 + A 2 [(dx + (3 l dt) 2 + a 2 g AB (dy A + [3 A dt)(dy B + [5 B dt)} (5) 

Here we use coordinates (x, y 2 , y 3 , t), where x and y A are spatial periodic coordinates on 
T 3 , with d/dy A being Killing fields, and t is a time coordinate (possibly, but not necessarily 
CMC). The functions a, A, /3 1 and (3 A are all functions of x and t (periodic in x), as is 
the two-dimensional unit determinant Riemannian metric cjab'-, the function a depends on 
t only. We also require that the density function \i depend only on x and t, and that the 
matter velocity field u a take the form u = vA~ 1 d/dx + wd/dt, with v a function of x and 
t only and with w determined by the condition g^u^u^ = —1. 

Substituting the metric (5) into the Einstein- matter equations (la), (2) and (3), and 
using the assumptions just stated, one obtains an explicit formulation of the Einstein-dust 
Cauchy problem for a class of T 2 -symmetric cosmologies in terms of a, A, /3 1 , (3 A , a, fj, and 
v. While we will need to use the local existence and stability results for the full Cauchy 
problem - see Proposition 1 above - we shall only be working with a couple of the equations 
in explicit form. The two we will need - they are constraints on the choice of initial data 
- are [12]: 

d 2 /dx 2 (A 1 / 2 ) = -1A 5 / 2 H(K - \t) 2 - It 2 + a 2 + 16^(1 + v 2 )} (6a) 
dK/dx + 3A~ 1 (dA/dx)K - lA- 3 d/dx(A 3 r) - A = 8nA/j(l + v 2 ) 1/2 v, (66) 

where r is the mean curvature of the chosen initial data (generally a function of a, ft 1 , fi A , 
a and their time derivatives) , where K is an eigenvalue of the second fundamental form of 
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the initial data - specifically one has: 



K = la' 1 [a~ l da/dt + dfi/dx] + \t 



(7) 



and where a and A are certain functions of the initial data which we need not specify, since 
they vanish for the cases we consider below. 

3. Initial Data for Spacetimes with Finite CMC Slicing 

In our examples the mechanism which is envisaged as stopping the CMC slicing from 
proceeding is the formation of 'shell-crossing singularities'. It is appropriate at this point 
to discuss this feature of dust spacetimes in some more detail. The general idea is that self- 
gravitating dust particles have a tendency to collide with each other, creating singularities. 
More specifically, in a spacetime with an isometry group with two dimensional orbits, a 
shell of dust particles which are related to each other by the symmetry moves in a coherent 
way. If two of these shells collide then the intermediate shells are trapped between them, 
so that the matter density is forced to blow up. Let us give a formal definition in the case 
of T 2 symmetry which is of interest here. 

Definition. Consider a T 2 -symmetric cosmological spacetime which is a solution of 
the Einstein- dust equations on which is defined a time coordinate t with compact level 
surfaces whose range is (£_,£+). Suppose that there is a time to with t- < to < t+ and 
two distinct dust particles which move orthogonal to the group orbits with the property that 
the distance between their world lines, as measured in the hyper surfaces t= const, tends to 
zero as t — > t + . Then we say that the spacetime develops a shell- crossing singularity as 
t -> t + . 

In the examples which we discuss in the following we will not prove that a shell-crossing 
singularity develops. However the strategy of the proof is guided by the idea that that is 
what is happening. To construct these examples, we rely on the initial value formulation, so 
our aim is to find simple initial data sets for Einstein-dust cosmological spacetimes, which 
are apt to lead to shell-crossing in the future, and which of course satisfy the Einstein 
constraint equations. We do this using the conformal method [5], as adapted to the T 2 
symmetric metric (5) and the corresponding matter variables discussed in Sect. 2. As in 
[5], we proceed as follows: 

- Choose (for some fixed constant to < 0) 




ii) 9ab = $ab (so A = 0) 

iii) a = 1 

iv) second fundamental form so that a = 

v) p,(x) any smooth positive function on 5" 1 with 




(8a) 



v(tt — x) = —v(x) 



(86) 
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- Determine K(x) by solving 

dK/dx = 8717/(1 + v 2 )v (9a) 

and A(x) by solving 

d 2 /dx 2 (A 1/2 ) = -\A^ 2 [lA-\k - \t ) 2 - \t 2 + 16nA~^(l + v 2 )} (96) 

- Set 

H = A~ 4 fl (10a) 
v = v (106) 
K = \t + A- l {K- \to) (10c) 

For a given choice of to, P>(x) and this method produces a solution of the constraints 
so long as equations (9a) and (9b) admit a solution K and A. The reflection symmetry 
conditions (8a)- (8b) guarantee that the right hand side of (9a) satisfies 

/i(l + v 2 )vdx = 0, 

s 1 

from which it follows that (9a) admits a solution (unique up to a constant). Sub and 
supersolution techniques [8] show that (9b) admits a solution as well, with that solution 
bounded below and above by the subsolution A_ and the supersolution A + respectively: 

A_ =t- 1 (24nB) 1 / 2 (11a) 
A + = Max{t - 1 (487r||A|U(l + ||i5 1|^) 1/2 , t~ 2/3 (|||K - ^o|lL) 2/3 } (H6) 

Here B := Min p, > 0. 

Note that, a priori, (3 1 and a may be chosen freely. If one wishes to maintain a CMC 
foliation into the future and past, however, then a must be chosen to satisfy 

d 2 a/dx 2 + A-\dA/dx){da/dx) = aA 2 [^(K - ±t ) 2 + \t 2 Q + 4^(1 + v 2 )] - A 2 (12) 

Similarly, to maintain a fixed x coordinate range, say [0, 27r], one must restrict (3 . We 
shall presume that this has been done. 

So far, we have not encoded shell-crossing into the initial data. This is readily done: 
one chooses a pair of points x\ and X2 in the interval (0, tt), and one chooses the function 
v so that v(xi) = 1 and v(x2) = —1- These data describe a pair of dust particles (among 
others) starting at x\ and X2 at to? with initial spatial velocities +1 and —1. If one uses 
CMC time - r(t) = t - then the equations of motion for these particles are 

dx/dt = aA- 1 v(l + v 2 )- 1/2 -P 1 (13a) 
dv/dt = -A^da/dxil + v 2 )~ 1/2 + aKv (136) 
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From these equations, as well as from the Einstein-dust equations, one obtains upper and 
lower bounds for changes in the velocities of the particles in finite CMC time [12]. Since 
these bounds are independent of \x2 — xi\, one expects that for sufficiently small \x2 — xi\, 
shell crossings are inevitable. 

Now, as a consequence of local existence theorems for the Einstein-dust PDE system, 
we know that for every choice of sufficiently smooth initial data, a cosmological spacetime 
consistent with those data - i.e., a cosmological development of that data - exists for at 
least a finite proper time into the future and into the past of the initial Cauchy surface at 
to- In the special case of T 2 -symmetric spacetimes, local existence theorems in terms of 
CMC coordinates have been proved [11]; so, one knows in addition that there is, in every 
maximal development of initial data constructed as above according to conditions (i)-(vi) 
and equations (9)-(10), a finite CMC time slicing in a neighbourhood of the initial surface. 

While one expects shell-crossing to stop a cosmological development from proceeding 
in proper time, it is possible a priori that CMC time could proceed to its limits - which 
are and — oo in a cosmological spacetime with E 3 = T 3 - in such a spacetime. As shown 
in [12], this is not the case. Indeed one has 

Proposition 2. [12] Let to < and e > 0. There exist sets of CMC (r = to) initial data 
for the Einstein-dust system such that in any development of the data, the CMC slicing 
which exists in a neighbourhood of the initial surface cannot be extended past either t + e 
or t — e. 

The basic idea of the proof is to use the Einstein-dust equations - including the CMC 
equation (12) - to obtain controls on the geometric quantities a, A, and K in a finite CMC 
time interval about to, and then argue that for \x2 — xi\ small enough, the dust particles 
starting at x\ and X2 (with initial velocities +1 and —1 as described above, and with metric 
described by equation (13)) must intersect within CMC time to + e, if the solution exists 
that long. An intersection of this type would violate regularity and so CMC time must 
stop. Similarly, if one chooses v(x) with v(xi) = —1 and v{x2) = +1, the intersection 
occurs in the past, within t — e. For details, see [12]. 

We note two facts regarding this result which are important for our work in the next 
section. First, we note that shell-crossing is essentially a local phenomenon. Specifically, it 
follows from the details of the proof of Proposition 2 that regardless of what is happening 
outside of the interval [xi,^], one can choose data inside [^1,^2] which results in shell- 
crossing occurring arbitrarily soon. Second, it is a trivial consequence of Proposition 2 as 
stated that given any sequence of positive numbers T m converging to zero, one can choose 
a sequence of Einstein-dust initial data C m such that in any spacetime development of C m , 
CMC slicing extends no further (in CMC time) than to + T m . 

4. Spacetimes Extending Past the CMC Slicing 

We now consider sets of initial data of the form described above, with shell-crossing 
built in, but with a 'quiet region' far away from [xi,^]- In this quiet region, we 

choose 

p,(x) = 1, X- < x < x+ (14a) 
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v(x) = 0, x- < x < x + (146) 

It follows from (9a) that (dK/dx)(x) = for X- < x < x+, so making a convenient 
choice of constant, we have K(x) = ^to on that interval. Thus we find that, modulo the 
conformal factor A 2 , which solves (9b) on S 1 , we have in initial data for a spatially 

flat Friedmann-Robertson- Walker cosmology. 

Our claim is that for all sets of initial data of this type, regardless of how soon shell- 
crossing occurs, A is sufficiently well-behaved in that in the maximal spacetime 
development of a given set of such data, the domain of dependence of extends 
a finite proper time into the future and into the past, and hence extends past the CMC 
slicing, if that is cut off sufficiently quickly. This is the content of our main result: 

Theorem 3. There exist Einstein-dust cosmological spacetimes, maximal developments 
of Cauchy data, which contain a CMC Cauchy surface but cannot be foliated by a CMC 
slicing. 

Proof. We fix an initial time to G (— oo, 0), and on the circle we fix a pair of disjoint 
intervals* [a?i,X2] (the shell-crossing region) and (the quiet region). We restrict 

attention to sets of data which satisfy conditions (i)-(vi) in the above procedure for con- 
structing initial data and equations (9)-(10), but with conditions (8a) and (8b) on p,(x) 
and v(x) replaced by the conditions that: 

a) p,(x) is a positive smooth function on S 1 with 

jj(x) = 1, X- < x < x + (15a) 

and 

fx(iT — x) = p,(x) (156) 

b) v(x) is any smooth function on S 1 with 

|v(x)|<l (16a) 

v(x) =0, X- < x < x + (166) 

v(n — x) = —v(x) (16c) 

v(xi) = 1 (16d) 

and 

v(x) = -1 (16e) 

for some x G (xi,X2). As discussed above, we can solve (9a) for K, and we obtain that 
K = |£ on [x-, x + ]. 

Now in accord with the comments at the end of Sect. 3, it follows from the details 
of the proof of Proposition 2 in [12] that one can choose a sequence of initial data C m of 
the form just described - with p,(x) the same for all m - such that the CMC slicing in any 



* Without loss of generality, we assume that [xi, x^\ fl [ir — x + , n — x_] = <f> as well as 

[x 1 ,x 2 ] n [x-,x + ] = 4> 
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development of C m extends no further than to ± T m . What we now wish to show is that 
A and its derivatives restricted to [x_,x+] are controlled uniformly, for all sets of data in 
this sequence C m . 

First, we recall the sub and super solutions (11) for A on S . The sub solution is 
clearly the same for all sets C m . For the super solution A^. : while ||/i||oo <md ||t>||oo are 
independent of m, the quantity \\K — |to||oo ma y n °t be. However, it follows from (9a) 
that 

\K - \t \ = 8tt| / u(l + v 2 )v\ 

3 ly^i (17) 

<i6 7 r 2 ||/i|| 00 (i + HiiiiDNlL 

which is independent of m. Hence we have upper and lower positive bounds for A(x) for 
the entire sequence. Next, as detailed in [12], we use a generalization of an estimate of 
Malec and O Murchadha [10] to argue that dA/dx is also uniformly bounded along the 
sequence. 

We now note the form of the Lichnerowicz equation (9b), which expresses d 2 A/dx 2 
in terms of A, /x, K, v and r. We know from above that A is uniformly bounded, and 
we have chosen jl and r to be fixed (and bounded) independent of the sequence. The 
only functions that change as m grows, and could cause trouble for d 2 A/dx 2 are v and 
K. However, in [x_, x+] we have v = and K — |to f° r all m; it follows that in [x_, x + ] 
we have d 2 A/dx 2 uniformly bounded along the sequence. Taking spatial derivatives of 
the Lichnerowicz equation (9b), we may similarly argue that spatial derivatives of A of all 
orders are, at least in [x_,x + ], uniformly bounded. 

With these controls on A, we have established uniform C°° bounds on the sequence C m 
of initial data, restricted to [x_,x+]. We may now apply the Cauchy stability statement 
of Proposition 1, from which we infer the following: there exists a pair of finite positive 
numbers A and fx (with A < ^(x + + x-)) and there exists a fixed harmonic coordinate 
system (x, y 1 , y 2 , t) on the spacetime region 

M 4 = (x- + A, x+ - A) x S 1 x S 1 x (t -u,t + v) 

such that for each set of initial data C m , there is a spacetime (M 4 , g m ) with the properties 

(i) (M 4 , g m ) is a solution of the Einstein-dust equations 

(ii) the metric g m induces the initial data C m on the hypersurface t = t 

(iii) the metric coefficients of g m and g^ 1 written in terms of the coordinates (x, y 1 , y 2 , t) 
are all C°° uniformly bounded 

In other words, because of the 'quietness' built into each set of initial data C m on 
the region [x_,x+], the maximal spacetime development g m of each C m lasts at least for 
a common (harmonic) time v. Note that we use harmonic coordinates here because the 
Cauchy stability proposition is proven [4] using harmonic coordinates. 

Let us now fix a point (x*, yf) E (x_ + A, x + + A) x 5" 1 x S 1 ; for each metric g m , we 
consider the timelike geodesic 7 m which is orthogonal to the initial surface 

(x_ + X,x+ + X) x S 1 x S 1 x {t } 

at the point (x*, y^ 4 , to)- Since g m and its inverse are both well-behaved in harmonic 
coordinates, one readily determines that for each j m , there is a positive number L m which 
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measures the proper time length of 7 m from the point y^, to) on the initial surface to 
the point where 7 m leaves M 4 . Moreover, the uniform boundedness of g m and g m x implies 
the existence of some L > such that L m > L for all m. 

By construction, in any spacetime development (M 4 , g m ) of the initial data C m , there 
is a CMC slicing in the neighbourhood of the initial surface, with the slicing lasting no 
longer than T m in CMC time and with lim™,^^ T m = 0. If we consider the length £ m of 
7 m restricted to the CMC region, we obtain the formula 

Cm = ! " + * a{n(t))dt (18) 

where t is CMC time. 

From the CMC equation (12) for a (which holds for all CMC t E [t ,t + T m ]) we 
note that at any point where a achieves a maximum, we have 

> A 2 a[\(K - \tf + \t 2 + 4tt//(1 + v 2 )] - A 2 (19) 

which implies (recalling that t < 0) that 

a < Z/tl (20) 

Hence from (18), we have 

Cm < (S/t 2 )T m (21) 

But linv^oo T m = 0, so that lim m ^ 00 C m = 0. This tells us that for some sufficiently 
large value of m, we have in M 4 (and in any development M 4 containing M 4 ) C m < L m . 
Thus in M 4 , the geodesic 7 m leaves the region of CMC slices before it passes out of the 
region with harmonic coordinates. 

Note that the arguments of this section do not rely essentially on the fact that the 
matter model is dust; a perfect fluid with pressure would do just as well. The only point 
where the fact that the matter model is dust is used is in applying Proposition 2. If 
an analogue of Proposition 2 could be proved for a perfect fluid with pressure then a 
generalization of Theorem 3 to that matter model would follow. It seems plausible that 
an analogue of Proposition 2 does hold for a fluid with pressure, due to the formation of 
shock waves, but the proof would be very different from that in the case of dust. On the 
other hand, there are other forms of matter for which the analogue of Proposition 2 is 
false; this has been proved for collisionless matter and wave maps in [11] 
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